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ABSTRACT 


The report presents successful results of determining total solar 
energy reaching the ground over the United States using measurements 
from meteorological satellites as the basic data set. The methods of 
satellite data processing are described. They are based on many years of 
prior experience, as well as harctvare and computer programs, related to 
quantitative analysis of meteorological satellite data. Uncertainty 
analysis and comparison of results with well -calibrated surface pyranometers 
are used to estimate the probable error in the satellite-based determina- 
tion of ground insolation. It is 10-15S5 for daily information, and 
about 5% for monthly values. However, the natural space and time var- 
iability of insolation is much greater than the uncertainty in the 
method. The most important aspect of the satellite- based technique is 
the ability to determine the solar energy reaching the ground oyer small 
areas where no other measurements are available. Thus, it complements 
the widely-spaced solar radiation measurement network of ground stations. 
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1.0 Introduction 


This Final Report under Grant NAS5-22372 presents results of 
research during the sixteen-month duration of the grant: June, 1976 - 
October, 1977. Several earlier reports were submitted as required. 

Support for the research was suddenly terminated when NASA's Office of 
Fnergy Programs was excised from the fiscal year 1977 budget. Despite 
this fact, the tedious development of the satellite data processing methods 
and computational algorithms have demonstrated feasibility of the satellite 
based method. Interest by both scientists and solar energy applications 
specialists has been extraordinarily high as measured by requests for 
the preliminary results of the study and by letters and calls inquiring 
about the method. In our opinion the positive results presented in 
this report indicate the advisability of continued research and demonstra- 
tion of methods for determining the solar energy reaching the ground 
using satellite data. 

Atmospheric scientists have studied solar energy for many years 
because of its interaction with clouds and gases in the atmosphere, and 
because of the surface heating due, in part, to solar energy reaching the 
ground. Both effects influence the thermal and dynamic state of the 
atmosphere and are important for weather forecasting. Indeed, the U.S. 
network of solar energy measurements was started and is maintained by 
the National Weather Service at their stations. 

With the advent of meteorological satellites in the 1960's, and 
the establishment of NOAA's operational satellite program shortly 
thereafter, atmospheric scientists used these new platforms to obtain 
measurements of solar energy reflected from the earth and atmosphere. 

The initial interest focused on the atmospheric problems noted above 
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(Hanson, et al,, 1967i Vonden Haar and Hanson, 1969). Throughout the 
course of this work the radiation specialists working on solar energy 
applications anci on atmospheric problems kept in close contact through 
such groups as the Committee on Radiation Energy of the American Meteorolog 
ical Society. The first author is a former chairman of this committee. 

At a special meeting to foster such coordination, sponsored by NSF (RANN) 
and organized by NOAA - the Solar Energy Data Workshop, November, 1973 - 
we proposed the method to determine the solar energy reaching the ground 
using measurements from meteorological satellites (Vonder Haar, 1974), 
During the following year, the research method was discussed and refined 
by the authors and their colleagues at the Department of Atmospheric 
Science, Colorado State University. At CSU, in the same College of 
Engineering, a larger, parallel effort by engineers including George 
Lbf and Susumu Karaki developed a nationally recognized program of 
solar energy applications research. We continue to benefit from dis- 
cussions with these engineers and scientists. Initially, the research 
leading to our larger effort with NASA was supported by a small grant 
from the State of Colorado, through its Colorado Energy Research Institute. 

During the 1974-75 period the solar energy data needs of users were 
specified in increasing detail. Many groups participated and the Report 
and Recommendations of the above mentioned Solar Energy Workshop (Turner, 
ed., 1974) contains typical results. In general, information is required 
about the (X,Y,T) variation of both direct solar radiation and total 
(direct + diffuse) radiation. The concentrators require input about the 
direct, while collectors (e.g. used for heating and cooling systems) need 
total solar information. Certain spectral regions are of special interest 
for photovoltaics, yet most information demanded by users can be supplied 
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by total spectral measurements for close approximation. The statistics 
of the energy reaching the ground are often as important to users as 
the mean value data. Variability of solar energy over the U.S. has been 
studied using the past records collected at the few weather stations having 
good continuous measurements from either pyranometers (radiometers) or 
sunshine recorders (see for example Pack and Korshover, 1974; Angel 1 and 
Korshover, 1975). The solar energy supply is characterized by both high 
spatial and temporal variability. Time scales from hours to years show 
variability. 

It is precisely because of this high variability that satellites 
are important measurement platforms for assessing our solar energy 
resource. Some important considerations are; 

1) Clouds are the principal modulators of the solar energy reaching 
the ground; meteorological satellites are designed to measure 
clouds. (See Figure 1.1 for a typical satellite photograph over 
the U.S.) 

2) The satellite radiometers have high ground resolution, from 1 
to 6 km, and are thus capable of resolving small regions of 
special interest (see Figure 2 in Appendix E for a high resolution 
view of the Lake Superior/Duluth area). 

3) Output from the satellite radiometer is digital (quantitative), 
thus amenable for computer processing (see again Figure 2, 

Appendix E). 

4) The satellite data are available from operational satellites, 
an integral part of NOAA's weather observation system, as fre- 
quently as every 10 to 30 minutes. 
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5) The entiire U.S. is viewed with the same Instruinent reducing 
problems of Intercalibration of many Instruments (e.g. at 
ground stations) and thus facilitating the Intercomparison of the 
Insolation between different sites, 

However, It Is a major step between the availability of the satellite 
Images and their potential yield of solar energy Information, The data 
processing system and algorithms noted In the following sections and the 
appendices to this report were the results of a lengthly research effort. 
This processing system for determination of solar energy reaching the 
ground considers several special aspects of the satellite data. Calibra- 
tion and location of satellite data were based on more than ten years 
experience In such research. The Important variable of surface reflectance 
and attenuation of solar energy by dust and water vapor In the atmosphere 
had to be Incorporated Into the processing system. Here we’ve drawn upon 
experience In the use of the weather data and on special techniques and 
computer programs to account for radiative energy transfer through the 
atmosphere. Finally, the sheer volume of satellite data dictates that 
special hardware methods for processing data be used. Most of these were 
also available from previous meteorological satellite research at Colorado 
State University, and were adapted for use In the solar research project. 

Such a major undertaking was prefaced by both possiblllty-of-success 
and worth-of-success considerations. Our initial proposal to NASA, based 
upon work sponsored by the State of Colorado, included a complete un- 
certainty analysis of the proposed method. We estimated an uncertainty 
of less than 9-10% for the determination of the magnitude of solar energy 
reaching the ground. Of course, we tould be sure that the spatial varia- 
tion, or gradient Information about insolation would be of very high 
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accuracy even if the absolute values were offset. (Results of our 
subsequent research* as described In the following sections, showed this 
to be a good estimate.) The Solar Energy Data Workshop Report contains 
a recommendation that solar radiation reaching the ground be determined 
with less than ± 5 % systematic error. However, the discussion of this 
goal at the workshop was extensive and Included the points that engineers 
are presently using data, where available, containing 10-15* systematic 
error. From the point of view of the satellite-method, our error estimates 
were not based on the use of corrections to the satellite- Inferred values 
by a high quality set of measurements at the surface . This Is certainly 
part of the proposed concept and thus we should do even better than the 
estimates would Indicate. Namely, the satellite-determined output will 
ultimately be adjusted to the new high quality data presently available 
from N0AA*s "baseline" stations at Mauna Loa, Samoa, and Barrow (Hanson, 
1971) and to the new, upgraded network over the United States (Riches, 
et a1., 1975). In summary, possibllity-for-success of the proposed 
method is found to be high. 

What about worth-of-success? At present and In the foreseeable 
future we have only about 35 to 40 high-quality solar energy measure- 
ment sites over the United States. When all lower accuracy radiometers 
and sunshine recording sites are Included, we have at best 100 data points 
across the entire country. Many of these are so Influenced by local con- 
ditions (e.g. urban effects) that they are unrepresentative of other 
locations only a few miles away. In addition. Pack and Korshover (1974) 
report that the spatial variation of monthly averaged sunshine between 
pairs of sites (e.g. Boston-Albany, Kansas City-St. Louis, Los Angeles- 
San Diego) can reach values of 20-30^4. Of course, daily spatial variation 
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can be 300% because of weather conditions. Thust how can the design 
engineer for a solar application project use such data frer?5 our surface 
network to consider a construction at a location without measurements? 

He Interpolates only with great uncertainty or more prudently, designs 
an excessive margin into his solar energy system. Some studies have 
estimated the cost of such excessive margins. Houghton (1976) considered 
three sites and modeled the economic impact of uncertainty in insolation 
data. Using monthly average statistics from the past data set he deter- 
mined that the sensitivity of annual fuel savings was about one dollar 
per year per percent difference in insolation multiplied by the homeowner's 
fuel cost per MBTU. Thus, better information, as from the satellite method, 
that reduces insolation uncertainty by 10% would save each homeowner ap- 
proximately $55 per year for a fuel cost of $5.50 per MBTU. In a separate 
study the Solar Energy Data Workshop Report states that users estimated 
benefit- to-cost ratios of 700 to 800 per unit cost of new ground station 
radiometers. Better information about the average solar energy supply 
and more extensive data on variability are obviously of substantial worth. 

Is the option of a tremendously increased ground-based network a 
real alternative? Instrument manufacturers have responded to the solar 
energy measurement need by discussing such proposals (Klein and Hickey, 

1973; Katz, 1975) and manufacturing a variety of new solar energy radio- 
meters. Some of these are of great importance to upgrading the few 
existing stations. However, the cheapest system which compares well 
with more expensive instruments (Silva Dias, et al, 1976) still entails 
about $500 per site for instrument plus digital recorder. Site acquisi- 
tion, maintenance and calibration costs, all labor intensive, are additional. 
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In our opinion* the optimum approach to obtain the necessary 
and valuable solar radlat-^cn microclimate over the United States, and 
to keep updating the data set. Is to: 

a) maintain a limited number of very high quality surface 
radiation measurement sites; 

b) fill the spatial gaps between such stations with a determination of 
solar energy reaching the ground based on satellite measurements. 

The following sections describe the satellite method, the data sets 
already available, and results that demonstrate feasibility of the method. 

2.0 Method of Determining Solar Energy Reaching the Ground Using Satellite 

Data 

After Vender Haar (1974) outlined the proposed methods, we refined 
the techniques to determine both total and direct beam solar energy 
reaching the ground In a set of papers (Vender Haar and Ellis, 1975, 

Ellis and Vender Haar, 1976). These papers are reproduced In their 
entirety In Appendix A to this report. As noted In Section 1.0 the 
methods are the results of many years of research using satellite 
measurements an ' radiation calculations to solve atmospheric research 
problems. Their adaptation to the determination of solar energy reaching 
the ground was a major result of the research presented In this report. 

3.0 The Satellite and Weather Data Used In the Research 

A principal objective of this study was to demonstrate a capability 
for determining total solar flux (insolation) at the ground, using earth 
orbiting satellite measurements of upwelling solar radiation flux. Thus, 
emphasis Is placed on applying satellite data to the problem rather than 
depending largely on a broad based ground measurement *hetwork. 
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Satellite data which qualified in various degrees of "goodness" tor 
the study are shown in Table 3.1. Some data sets are archived strictly 
as photo images and others as digital radiance values. All satellites 
listed in the table were in sun synchronous orbits sampling at just one 
local daylight sun time each day with the exception of SMS 1 and 2, An 
SMS satellite has a geostationary orbit such that the satellite subpoint 
is nearly fixed at one latitude and longitude. Each area element is 
sampled approximately every 1/2 hour. Most of these data are in imagery 
format with only a few intermittent digital sets being archived. 

Satellite data actually used in this study are as follows: approximately 
120 days of ESSA-^9 imagery concurrently with Nimbus 3 digital radiance dai-u 
for May 1969 through January 1970; daily digital data for the month of 
June 1975 from NOAA-4. (Some digital and Imagery data from SMS 1 and 2 
were acquired but were not processed because of a premature discontinuance 
of grant funding.) 

To establish an empirical function for atmospheric absorption in 
cloudy atmospheres we required daily U.S. rawinsonde network observations 
and a reliable ground network of daily pyranometer measurements. For 
the June 1975 period there were just 9 sites in the National Weather 
Service ground network which recorded total daily Insolation with well 
maintained Eppley Model II pyranometers. Absolute accuracy of this 
pyranometer is +, B%. The 9 stations are shown in Table 3.2. 

Daily ROAB's are also required so that computation of daily precipitaDle 
water (W) can be made. The W values are computed at each RAOB site and 
then objectively analyzed at 1000, 850 and 700 mb levels to fit the NWS 
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TABLE 3.1 
SATELLITE DATA 


SATELLITE 

TIME 

GRIDDED 

SPACE 

IMAGERY' 

DIGITAL^ 

EssA y 

1430 

6 NM 

* 


NIMBUS 3 

1130 

200 KM 


•k 

ITOS 1 

1430 

6 NM 

* 


NOAA 2 

0900 

6 NM 


k 

NOAA 3 

0900 

6 NM 

★ 

k 

NOAA 4 

0900 

6 NM 

ie 

k 

SMS 1 

1/2 HR 

1/2-2 NM 

★ 


SMS 2 

i/2 HR 

1/2-2 NM 

★ 



Vhe * symbol indicates the format of the archival (Circa, 1976). 
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TABLE 3.2 

STATIONS WITH EPPLY MODEL II PYKANOMETER DATA FOR JUNE 1975 

Daily Hourly 

Albuquerque, NM * * 

Bisiaarck, ND * 

Caribou, ME *• * 

Ely, NY * 

Las Vegas, NV * 

Miami, FL * * 

Oak Ridge, TN 
Phoenix, AZ 

Seattl e-Tacoma, WA * 
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forecast grid spacing Which Is every 32nd fine mesh point (approKlinately 
10 km on a side) of the satellite data grid. The generated M fields 
are» at a later step In the final Insolation computation, Interpolated 
in three dimensions to the fine mesh satellite data grid. A complete 
description of the various operations performed on each data set and how 
the data sets come together for a surface Insolation computation appear 
in Appendix B. 

4.0 Results of the Study 

In this section the atmospheric absorption parameterization is 
discussed along with derived values of surface Insolation from the satel- 
lite measurements. Following that, we present a comparison of satellite 
results with Insolation climatologies from a ground network of measurements 
for a test of the satellite method. 

4.1 Absorption Parameterization 

As discussed in Ellis and Vender Haar (1976) (in Appendix A-2) the 
fractional atmospheric absorption of solar energy (q^) in a cloudy at- 
mosphere was found by Hanson (1971) to have an empirical relationship 
to optical pathlength (t) in the form 

pg « Cl + Cg Inx (4.1) 

Hanson (1971) derived coefficients Ci » 0.117 and Cg ® 0.31 by 
applying a least squares fit to q^ and t, both of which were determined 
at a ground pyranometer site. Satellite data used in his study consisted 
of monthly mean values of mediuL ii;:solution (50 km at satellite subpoint) 
Nimbus 2 scanning radiometer measurements. The value x is a function of 

^In this method the clouds scatter, but do not absorb the solar energy. 


■ s* X 
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amount of cloud (measured from satellites) and water vapor optical depth 
(from RAOBS). 

The daily and t values (determined from the method presented in 
Appendix B) for the nine ground based pvs'ao»'ometer locations for the 
month of June 1975 are plotted in Figure 4.’!. The ground based set was 
limited in size to include only the well calibrated and maintained 
instruments within the NWS network of global Insolation measuring stations. 
The curve in Fig. 4.1 is a least squares fit to Equation 4.1, the 
coefficients being Cl = 0.133 and C2 = 0.033. 

There is a considerable amount of scatter of the points about the 
line of regression. The scatter of the values at each of the nine 
stations about the nine-station regression line was computed and is 
shown in Table 4.1. The scatter of daily values Is on the order of 10 
percent of the insolation received at the top of the atmosphere. The 
nine-station parameterization of q^ is within two percent zt six of the 
stations and within six percent at the remaining three stations (ave*^age 
algebraic difference is defined as the q^ computed from Equation 4.1 
minus the computed at each station from Equation 1 in Appendix B 
averaged over the month). 

The Miami and Seattle data sets were extracted from the nine station 
data set in order to assess the scatter of points observed in Fig, 4.1. 

The ^ stations were selected as extreme cases of departure of q^ from 
the line of regression, i.e., the scatter at Miami is 13.6 percent and 
the average algebraic difference at Seattle is 6.0 percent. The plot of 
q, versus t for the two stations is shown in Fig. 4.2. A least squares 

a 

line of regression of the form of Equation 4.1 was fitted to each data 
set in order to aid comparisons. The two curves are similar but differ 
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Figure 4. 1 


Least squares fit for daily values 


over 9 pyranometer 


sites. 


ORIGINAL PAG€ IS 
OF. POOR QUALITY 
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TABLE 4.1 

Differ*ences between atmospheric absorption derived from a regression 
function and that derived as a residual of. measurements. 


Avg. Algebraic Diff Scatter 

Station Percent ly day“^ Percent ly day"^ 


Albuquerque 

-1.3 

-13 

6.4 

63 

Bismark 

0.5 

6 

8.4 

84 

Caribou 

0.4 

4 

8.1 

81 

Ely 

O 

1 

-40 

9.6 

96 

Las Vegas 

2.4 

24 

7.1 

70 

Miami 

-1,5 

-15 

13.6 

131 

Oak Ridge 

-1.9 

-19 

9.5 

94 

Pheonix 

0.9 

9 

6.0 

59 

Seattle 

6.0 

59 

9.7 

96 

Nine 





Stations 

0 

0 

1 

1 

CO 

CO 

42 


Figure 4. 



2 A daily plot of atinospheric absorption versus optical 
pathlength at two stations for the month of June 
1975. 
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considerably from the nine-station curve of Fig. 4.1. The mean vertically 

Integrated water vapor content of the atmosphere for June 1975 at Miami 

Is 4.2 cm while at Seattle it Is 1.6 cm. The larger precipitable water 

content at Miami accounts In general for the larger q. at Miami but does 

not account for the large slope and the offset between the curves of 

Fig. 4.2. The offset may be caused by diurnal changes In cloud amount, 

and thus q^, at Miami. The effect of a daytime average cloud cover (A^) 

greater than that observed by the satellite during its overflight would 

be to reduce q^ by increasing q^ and to reduce t (;;ee Appendix B). 

Further assessment of the diurnal cloud change effect requires knowledge 
3qr 

Appendix A2 shows that, q^, the surface Insolation away from the 
pyranometer location is determined from; 


‘’s = 


( 1.0 “ %-%) 

(1.0 - Pg) 


(4.2) 


We obtain q^ , as discussed above, and q^ and are measured from the 
satellite. 


4.2 Satellite Derived Surface Insolation Maps 

Prior to demonstrating a capability using the higher spatial 
resolution satellite data in June 1975 a study was done with medium 
resolution scanning data from the Nimbus 3 satellite (gridded at ap- 
proximately 200 km on a side area) for the state of Colorado only. The 
original coefficients of Hanson (1971) were used in equation 4.1 to 
determine atmospheric absorption. Figure 4.3 shows these results for 
six grid areas covering Colorado. It is quite obvious that this large 
grid spacing is not ideal because of the large variety of climate classifi- 
cations within each grid area, e.g. valley, mountain, plateau, desert, 
forest, and etc. 
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Results from the higher resolution solution NOAA 4 scanning radiometer 
data are shown in Figure 4.4. The gridded results of our study are of 
such high spatial resolution that they would not give a good basis for 
comparison to other studies. Thus» every sixteenth grid point in each 
dimension was selected for the isoline analysis in Figure 4.4, The 
general features are in qualitative agreement with what one would expect 
becaus^of cloudiness associated with various regions of the U.S. There 
are low values in the Pacific Northwest and New England States, and high 
values in Southwestern desert regions. (See Appendix E for further 
discussion.) 

4.3 Comparison With Other Insolation Climatologies 

The image in Figure 4.4 is compared with Figure 4.5, a mean June 
analysis derived from historical ground based measurements (U.S. Climatic 
Atlas). Both figures are in langleys per day. Since the latter figure 
represents an average of many months of June measurements while the 
former is but for June 1976, only qualitative comparison will be made. 

Agreement is seen in the two figures with low values in NW, NE and 
SE U.S. and with high values in West and Southwest U.S. A finger of low 
insolation extends northward into the Mississippi Valley which is re- 
stricted to the valley in the satellite results. Some disagreements are 
worth noting. A maximum appears in the satellite results over the 
Colorado Rocky Mountains which is undoubtedly a bias towards the morning 
hours at satellite observation time (9 a.m. local solar time) before 
midafternoon cloudiness begins -a frequent occurrence over the Rockies. 

The extension of the 600 ly/day contour in Figure 4.6 into NE Colorado 
was based on the Boulder, Colorado data - a station which is in the lee 
of the mountains and is subject to considerable mountain induced wave 
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Figure 4. From the Climatic Atlas of the United States, U.S. Department 
of Commerce, Environmental Science Services Administration, 
Environmental Data Service, June, 1968. 
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and orographic cloudiness. The coarse grid satellite results Indicate 
that the Isollnes to the east of Colorado are not to be extended west- 
ward Into Colorado. 

Mean dune Isollnes of Insolation, computed by Dennett (1965) 
applying multiple regression relationships between 12 years of Insola- 
tion and percent possible sunshine, cloud cover, and station elevation, 
are shown In Figure 4.6. Agreement with satellite results Is seen In 
the Great Plains and over Wisconsin. 

Examination of Figure 4.7, a compilation of Insolation statistics 
for a number of years at the ground by Texas A&M University, shows low 
insolation values in Eastern Texas which were not apparent In Figures 
4.5 and 4.6. The satellite »'asults show a trough In the 500 Isoline In 
Eastern Texas. The satellite results also show minima along the Gulf 
Coast of Texas. These minima are quite possibly due to anamolous 
cloudiness In June 1975. 

It Is quite unsatisfactory to compare maps of one month of June 
statistics (1975) to time averages of many months of June In any more 
detail than was done above. However we were able to assemble the data 
on "percent of possible sunshine" as actually reported by the National 
Weather Service for June 1975 (Figure 4.8). Comparison of their map 
with the totally Independent satellite results of Fig, 4,4 provides 
strong verification of the satellite method. Both the Minnesota and 
Gulf Coast minima detected by the satellite are shown In the surface 
sunshine data. Thus It does appear that except for regions which have a 
large daytime variability In cloud cover (Rocky Mountains, etc.), the 
sun synchronous, once-a-day satellite results are comparable to spatially 
smoothed results from point measurements of the grou>':d-based pyranometers. 


Figure 4.6 From “Monthly Maps of flean Dally Insolation for the United 
States" by I. Bennett. Presented at the Twelveth Conference 
on Radar Meteorology, Normati, Oklahoma, October 17-20, 1966 











PERCENTAGL OF POSSIBLE SUNSHINE, JUNE 1975 



Zooming in to use the satellite results at full spatial resolutions 
will give insolation statistics which would otherwise be available only 
frofii a spatially dense ground network of measurements (see example in 
Appendix E). The present study can be improved by eliminating errors 
due to daytime diurnal variations in cloud amount by applying the same 
methodology to quantitative SMS/60ES radiance data. These data are 
now available at the new Direct Readout Satellite Groundstation (Fig. 4.9) 
recently installed at Colorado State University. 
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Appendix A-1 


ABSTRACT 

A method is presented for determining solar insolation at the 
earth's surface using satellite broadband visible radiance and cloud 
imagery data, along with conventional "in situ" measurements. Conven- 
tional measurements are used to both tune satellite measurements and 
to develop empirical relationships between satellite observations and 
surface solar insolation. Cloudiness is the primary modulator of sun- 
shine. The satellite measurements as applied in this method consider 
cloudiness both explicitly and implicitly in determining surface solar 
insolation at space scales smaller than the conventional pyranometer 
network. 

INTRODUCTION 

Solar energy received at the earth's surface is viable as an addi- 
tional source of energy. There is now a need for more knowledge about 
its variability in space and time so the solar collector, concentrators, 
and storage devices may be designed for optimum performance. 
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Previously, surface measurements have been used to assemble regional 
to global spatial scale maps of monthly means of surface solar energy. 

These have been put together by various authors, some of which are: 

Fritz and McDonald (1949), Ashbel (1961), L5f, et al. (1966), and Bennett 
(1965). Cloudiness Is the primary modulator of solar energy. Satellites 
can observe cloudiness over large areas with more spatial Integrity than 
surface observations. Satellite measurements now available allow solar 
energy determination over a wide range of spatial scales much smaller 
than the conventional network permitted In the earlier works. Additionally, 
satellite measurements permit estimates to be made In areas where no con- 
ventional data are available. 

Broadband shortwave radiances (0.3 to 3pm) measured by the satellite 
along with cloudiness derived from satellite Imagery are used to compute 
surface solar energy. Conventional surface pyranometer measurements are 
used to "tune" the satellite measurements. 

SOLAR ENERGY AT THE EARTH'S SURFACE 

Solar energy reaching the surface of the earth beneath a 
heterogeneous cloudy atmosphere is a function of many variables. The 
irradlance at the surface, (energy per area per time) is the Integral 
over 2n steradians of the radiance reaching the surface, N (energy per 
area per time per solid angle). 

Ns(0,'t) = Q^cose^ e"'' ^ q„cos0„ f(e,P5) 

Energy reaching the surface may be either direct beam (first term) 
or diffuse "skylight" radiation (second term). Symbols in equation (1) 
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are defined as: 

" surface albedo 

■» solar Irradlance at "top of the atmosphere" (the solar 
constant) corrected for earth-sun distance 
0,<j) » zenith and azimuth angles, respectively 

t * optical depth of the atmosphere 
Q » solid angle subtendeo by the sun at the earth 

” function representing diffuse radiation as discussed by 
Diermendjian and Sekera (1954) 

A method for determining Q^, the total (direct-plus diffuse) radia- 
tion and a method for approximating the first term, the direct, follow. 


DIRECT PLUS DIFFUSE DETERMINATION 

The total broadband visible Irradlance received at the top of the 
earth's atmosphere, Is balanced by absorption and reflection within 
an earth-atmosphere column as In equation (2): 

Qq « RRy + RA|^ + RAg (2) 

where symbols are defined as: 

RRj = total reflected energy to space 

RA^ B total absorbed energy within the atmosphere 

RAg = total absorbed energy at the earth's surface 

These and other interactions of a single stream of radiation en- 
tering the atmosphere from space are shown In Figure 1. No secondary 
reflections are shown for simplification. 

An incoming radiation stream undergoes losses within the atmos- 
phere, other than absorption, due to reflections to space by clouds, 
RRq, by variable atmosphere particulates, RRp, and by background 
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levels of atmospheric constituents, such as by Rayleigh scattering, 

RRg. These three reflections make up the total atmospheric part of 
reflection to space, RR^. 

The downward directed stream at the earth's surface surviving 
atmospheric attenuations, is balanced by the sum of energy ab- 
sorbed at the earth's surface, RAg, and that reflected from the earth's 
surface, RRg. This relationship can be shown to be: 

RAg “ QsO-Ps) 

where is the earth's surface albedo, (p. * RR^/Qg). All terms 
shown in Figure 1 can be expressed as fractions of the Incoming beam 
Qq, and from now on will be shown as lower case letters, l.e., rr^ ** 
RRj/Qq, etc. 

Equation (2) Is rewritten using equation (3) and rearranged to 


give: 


% * 


1.0 - rrj “ ra^ 

0 - Pg) 


(4) 


Equation (4) is the relationship for fractional surface solar 
energy, q^, to which the satellite data will be applied. 

The term rr-j- is the actual planetary albedo at the top of the 
atmosphere measured by the satellite. The term pg, the surface albedo, 
cannot be measured directly from the satellite but can be estimated 
by using clear sky albedoes for the earth-atmosphere column. The 
clear sky albedo can bi? determined for each area, a grid area, by 
taking the minimum value measured at that area over some prescribed 
time period and number of samples. Implicit in the minimum albedo 
technique is that cloud effects (rr^,) and some variable atmosphere 
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particulate effects (rrp) are high frequency phenomena at a given grid 
area so that at least one sample over the prescribed time period Is 
free of these effects. The clear sky albedo, i^elated to 

the surface albedo, p^, by: 

<*‘’"T>in1n ■ ‘'s "s * ''V 

where is the transmittance of the atmosphere to a stream of energy 
reflected from the earth's surface and rr^f is the reflection to space 
by the atmosphere alone. 

Equation (5) can be solved for p^: 


D “ 
*^S 


Is 


( 6 ) 


The q. appearing in equation (6) is not the same as the q^ in 
the equation (4). All terms in equation (6) are for clear sky (cloud 
free) conditions. The terms q^, rr^t and tI can be estimated well 
from the nearest surface pyranometer station data with appropriate 
adjustment for any slight differences in solar zenith angle, 

Basically, then, the earth's surface albedo is determined from the 
clear sky albedo by "tuning" equation (6) with station data. 

The remaining term in equation (4) which must be satisfied is 
ra^, the cloudy atmosphere fractional absorption. Fortunately, ra^ 
is the smallest term in the expression (less than 0.25) since it 
cannot be solved for directly over the satellite measuring domain. 

It can be estimated sufficiently well by developing empirical re- 
lationships based on cloudiness and optical path length due to at- 
mospheric water vapor. This method is similar to Hansen, et al . (1967) 
and was later refined by Hanson (1971) by explicitly considering cloudiness. 
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He developed an empirical relationship for the month of June using 
conventional data so that; 

ra^ * 0*117 + ,031 t^ Ini (7) 

where t Is computed as a function of cloudiness and optical path length 

5 - 

I Of water vapor. The standard deviation of the residuals from least 

square fit In equation (7) was 3,2 percent. 

Small Scale Reduction . In addition to computing qg at the spatial 
scale of the satellite radiance data using equation (4), q^ can be 

y approximated at even smaller scales or at subgrid scales If cloud 

f amount 1s known at the smaller scale, Consider q» at the larger 

S a 

I scale so that: 

! q_ « a^ q^ + a_ q_ (8) 

I ^s c ^c n ^n ^ 

I 

I where a. and a„ are the fractional cloud and cloud-free areas over 

u n 

the larger scale, and q^ and q^ are the surface solar energy In the 
cloudy and cloud-free areas, respectively. Equation (8) can be solved 
for q^ provided the cloud amount 1s known, Essentially, q^ and q^ 
are average values of surface solar energy beneath an average cloud 
and in the clear regions, respectively. A measure of the distribu- 
tion of cloudiness over the small scale permits an estimate to be 
made of the distribution of surface solar energy over the area at a 
smaller space scale (subgrid scale). 

Cloud Area . Cloud cover as seen from satellites can be computed from 
satellite cloud Imagery. Figure 2 Is a photograph of cloud Images from 
the Defense Meteorological Satellite Program (DMSP) satellite. The 
cloud areas can be retrieved from this imagery by either digitizing 
the photograph or using an optical planimeter. These cloud areas should 
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be more properly termed as relative cloud amounts since the retrieval of 
true cloud amount is limited by the resolution of the satellite viewing 
instrument. 

The rectangularly enclosed region in the photograph in the left of 
Figure 2 has been digitized and is displayed to the right. The total 
digital signal has been grouped into five classes represented by just 
five different characters for display purpose. From the digital product 
relative cloud areas can be computed by selecting an appropriate digital 
cutoff value for cloudiness. 

This approach in getting subgrid scale solar energy at the earth's 
surface assumes that clouds over the entire large scale area have the 
same radiative properties. Quite often this will not be the case. 
However, this approach allows a best esy.1male of total surface solar 
energy at the smaller scale in lieu of direct ground measurements. 

Direct Beam Approximation . The direct beam solar energy is estimated 
with an assumption that clouds, other than thin transparent clouds, 
have a binary effect on It. That is, in the presence of a cloud, only 
diffuse solar energy reaches the earth's surface beneath a cloud. With 
this assumption direct reflections from adjacent clouds are not con- 
sidered. Thus, direct beam solar energy is received at the earth's 
surface only in the cloud free area. 

This method is to be applied over a small region, such as a city, 
using very high resolution satellite cloud imagery. For it to work 
well, individual clouds should be distinguishable in the imagery. A 
prior knowledge from surface measurements using pyranometer and pyr- 
heliometer instruments is required to separate the diffuse and direct 
solar energy in the clear areas. 
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If equation (8) is applied to the very high resolution data In 
regions where larger scale has been computed, then diffuse radiation 
at the surface In cloud areas can be estimated at the smaller scale too. 
Therefore, both direct and diffuse energies received et the earth's surface 
can be evaluated and statistics can be compiled on them. 

In the total and direct beam methods just described, the primary 
role of the satellite observation is to consider the effect of horizontal 
Inhomogeneities in cloudiness over an area rather than just scattered 
point values of insolation. Since clouds are the principle modulators 
of surface insolation, a solar energy climatology generated by these 
methods should be more representative than that generated by previous 
methods which were not augmented with satellite data. 
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INTRODUCTION 

The purposes of this paper is to demonstrate a method for determining 
mean daily solar insolation (direct and diffuse) received at the ground 
using satellite visible channel radiance measurements. The method Is 
basically an expansion of a method first demonstrated by Hanson (1971). 

An advantage of using satellite measurements is that one Instrument 
can take samples over a large region at fairly high spatial resolution. 

By using satellite data archives presently available, an insolation 
climatology can be generated even in remote regions at minimal cost. A 
principal disadvantage to the satellite method using archived sun- 
synchronous satellite measurements Is that just one measurement for each 
day at each ground position is taken, Thus, regions with significant 
diurnal cloudiness variations are not adequately sampled. However, 
using geosynchronous satellite measurements, such as those from a well 
calibrated SMS/GOES satellite, would eliminate the diurnal sampling 
problem. 
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A window, enclosing the contiguous United States, was extracted 
from the global archive of NOAA sun- synchronous satellite measurements 
for 30 days of June 1975 to demonstrate the method. The data are 
geometrically normalized Into albedos and combined with ground eleva- 
tion, precipitable water, and global pyranometer data to arrive at 
a high spatial resolution, mean dally Insolation for the month*. 


SATELLITE METHOD 


Total (direct-plus-diffuse) Incoming solar flux at the ground can 
be expressed In an energy balance equation of an earth-atmosphere 
column. The development of this expression Is taken from Hanson, et al. 
(1976). It Is normalized by dividing through by to get! 

1 = qr + <*a + qg 0) 

and 


q„ = (1.0 - pj q, 


( 2 ) 


where : 




Q /Q 
^x' ^0 

incident solar flux at top of atmosphere 
fractionally reflected shortwave flux 
fractional absorbed shortwave flux within the atmosphere 
fractional absorbed shortwave flux at the ground 
fractional incident solar flux at the ground 
ground surface albedo 


From (1) and (2) one can obtain such that: 
1.0 - q - a 


(1.0 - Pj) 


( 3 ) 
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Satellite visible channel radiance measurements are applied to 
(3) to obtain q^. Now is the actual planetary albedo which can 
be determined directly from satellite nteasurements with appropriate 
angular corrections. These corrections are discussed in the next 
section. The absorbed shortwave flux, q^, can be determined as a 
residual in (3) at sites where incident global solar flux is measured. 

Hanson (1971) discussed the development of an empirical formula at 
these sites which is then applied away from the sites. In effect, it 
allows one to extrapolate away from the sites using a statistical- 
empirical formula. His method is to compute values of total optical 
pathlength, t, and with determined q^ values, find coefficients to a 
functional relationship between the two using least squares. 

da “ + C 2 T “ Jin T (4) 

and considering diffuse radiation beneath clouds: 

•c = (T-Ac) (.86 U*+.14(1.66 u)) + A^(1.66 u) (5) 

where: 

A- = fractional cloud amount 

u = optical depth of water vapor for clear sky (precipitable water) 

u* = u/cosc,S is average, cosine weighted, daytime solar zenith angle. 

This basic method was expanded in the present application. The 
fractional cloud amount, A^,, is approximated from the satellite de- 
termined albedo rather than using scattered ground observations of 
A^. An obvious advantage to measuring A^ from a satellite platform 
is the consistent and very good spatial sampling capability of the satellite, 
even in remote regions of the earth. A method similar to that of Reynolds 
and Vender Haar (1976) is used to approximate A_. 


" %ld'%rr 
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( 6 ) 


where: 

‘^clr “ planetary albedo 

^cld “ Pl3inetary albedo over clouds (set at 0.50) 

They showed quite favorable results from many test cases for three 
different sites. 

The study was further expanded by approximating ground albedo Pg, with 
a cloud free albedo at the top of the atmosphere. In most cases the two 
differ because of atmospheric absorption and scattering. However, the 
difference is virtually never more than 5 percent in real atmospheres 
(Cox, 1975). This difference does not severely affect the results in 
the study as is shown by an uncertainty analysis in a later section. 
Furthermore, cloud free albedos, which are determined by retaining a 
minimum value of albedo at each grid spot on earth over some time period, 
are most likely more accurate than values extrapolated from rather 
sparse surface albedo measurements. 

DATA REQUIREMENTS 

Data needed to apply this method included calibrated visible- 
channel satellite radiance measurements, 12 hourly rawinsondes, and 
a number of well calibrated, global s51ar flux measurements at the 
ground. 

The method was tested on the month of June 1975 using NOAA 4 
scanning radiometer (SR) visible channel photometric measurements. 

The SR data from NOAA satellites are archived by the National Climatic 
Center as normalized brightness values. The archive extends back to 
1973 to include NOAA 2, 3, and 4 measurements. The measured brightness 
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values were normalized before archival by dividing each by cos ? (? is 
the solar zenith angle at time of measurement). Ground resolution of 
the archived data are about 10 km over the contiguous United States. 

The region where the method was applied is shown in a window ex-^ 
tracted for 11 June 1975 from the NOAA 4 global data set of normalized 
brightness (Fig. 1). Picture quality is reduced considerably because 
4x4 grid elements were taken together to obtain levels of gray shading 
on a microfilm plotter. 

The archived data must be converted to radiometric units from 
photometric units, and corrected for angular dependence of reflecting 
surfaces in the field of view. A model for correcting for angular 
dependence has been discussed in detail by Raschke, et al. (1973). It 
corrects for bidirectional reflectance characteristics of a grid element 
giving an integrated global flux over 2ir steradians. These instantaneous 
flux values are further normalized to mean daytime flux values by applying 
a reflectance model which allows for changing albedo with changing solar 
zenith angle. 

Angular corrections to be applied to the archived normalized brightness 
values in order to obtain mean daytime flux are shown in Fig. 2 for just 
one day, 1 June 1975. Three separate orbital swaths were mapped in this 
particular window. A discontinuity in correction factors is seen at the 
overlap boundaries along with an increase in values from left to right 
for each swath. 

Optical pathlengths due to water vapor, u, are computed from 12 
hourly rawinsonde observations within the window. The u values are 
objectively analyzed daily and interpolated to grid points, corresponding 
to the NMC prediction grid using a Cressman (1959) type scheme. The values 
at the NMC grid spacing are linearly interpolated to the NMC superset 
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grid spacing of the to the NOAA 4 archive. They are further adjusted 
for variations In ground elevation allowing for a linear decrease of u 
In the vertical with an Increase In ground elevation. Ground elevations, 
mapped to the nearest 10 meters on the MMC superset grid, are displayed 
In Fig, 3 for the window. 

Finally, Incident global solar flux at the earth's surface Is taken 
from the NWS .urface pyranometer network of measurements. From this 
network, sites which have well calibrated Epply Model II pyranometers 
were selected. Oust the 9 sites shown in Fig. 4 satisfied this criteron 
In June 1975. Remaining sites of the network were not included because 
of large uncertainty in their absolute calibration and relative measure- 
ment accuracy (see Hanson, 1973 for discussion of uncertainties in the 
network measurements). 

UNCERTAINTY IN THE RESULTS 

The method not only allows for high spatial global solar flux 
determination at the ground but also gives results at the ground com- 
parable in quality to in situ measurements. A complete uncertainty 
analysis of both independent and dependent measurement error shows an 
uncertainty in the results of less than 6 percent of the incident solar 
flux at the top of the atmosphere, or less than 8 percent of the clear 
sky value. However, these results assume no diurnal change in cloudiness 
and other atmospheric constituents. Of course, this assumption is not 
valid at many locations during various seasons of the year. Considering 
diurnal changes, an uncertainty closer to 10 percent in daily mean 
values seems reasonable. Since computations are made daily, statistics 
of various types other than just mean daily values can be generated. 

These statistics include variances and consecutive days of solar insolation 
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below some critical value. Furthermore, a number of additional products 
can be generated from the elements of the study. These Include a 
representation of cloudiness statistics and atmospheric absorption 
statistics. 

Finally, results of the study will be presented at the conference. 
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DATA ASSIMILATION AND COMPUTATION NETWORK OF EVENTS 
LEADING TO INSOLATION COMPUTATION FROM SATELLITE MEASUREMENTS 

Appendix D 

INTRODUCTION 

The solap energy microclimate study Involves quite a few different 
data sets along with a number of operations to be performed on each 
data set. It seems best to discuss the program by referring to a 
simple network of Identifiable events which take place for program 
completion. 

Three generalized networks are discussed in the following order: 
calibration and calibration verification network, data processing and 
mean Insolation network, and statistics generation network. Events of 
the networks are coded for reference to documentation. 

CALIBRATION NETWORK 

The NOAA satellite scanning radiometers are calibrated photo- 
metrically in luminance. The digital mapped data archive, mapped to a 
northern hemisphere polar sterlographic daily mosaic, is In the form 
of counts from 0 to 255. These counts are related linearly to foot- 
lamberts. 

Our research requires radiance data (Wm"^ sr"b* All calibration 
Information required to make the transformation from luminance to radiance 
has been received except for spectral radiance profiles of each calibration 

TV 

lamp. These profiles, representative of the lamp type used In the 
calibration, are most likely available. A representative from NOAA/NESS 

* A spectral radiance profile for the calibration lamp-system was received 
from A. Schwalb, NOAA (ref. Jones et al., 1965). 

Jones, G.D. , D.T. Hilleary, and B. Fridovich, 1965; A diffuse light 
source for calibrating meteorological satellite television cameras. 
Applied Optics . £, 3. 
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is trying to obtain this profile information for us from Santa Barbara 
Research Center. 

The original plan for both calibration and calibration verification 
is shown by the calibration network in Figure 1, Each event is discussed 
separately in numbered sequence as follows; 

1. Obtain calibration lar information. 

2. Obtain NOAA satellite scanning radiometer calibration test 
data sets, Four sets were obtained to compare with ATS-6 GVHRR 
visible channel data and one set to compare with Nimbus 6 ERB 
scanner data. 

3. Develop program to retrieve angles specifying geometry between 
satellite, sun and target. 

4. Apply calibration, angle retrieval program, and bi-directional 
reflectance model to NOAA-SR test set to obtain reflected flux 
and albedo. 

5. Obtain ATS-6 GVHRR data set to check NOAA-SR calibration. The 
ATS-6 data were not delivered from NASA/GSFC due to data re- 
duction problems. We proceeded without this calibration check. 

6. Apply bi-directional reflectance model to ATS-6 data (did not do). 

7. Compare NOAA-SR and ATS-6 simultaneous flux data for the selected 
target as a calibration check. One of the targets is White 
Sands Missile Range (did not do). 

8. Obtain Nimbus-6 ERB scanner data set to check NOAA-SR calibration. 

9. Compare NOAA-SR and Nimbus 6 flux data for selected targets as 
a calibration check on the NOAA-SR data. 

MEAN SOLAR INSOLATION NETWORK 

This network explains the sequence of events for computing mean 
monthly solar insolation at the ground. 
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The sequence of events and the status of each event is displayed 
in the mean insolation network of Figure 2, Each event is discussed as 
follows: 

1. Complete calibration, compute satellite-sun- target geometry, 
apply bi-directional reflectance model and evaluate reflected 
flux by comparison with other calibrated satellite flux data. 

This event is explained in detail in the previous section on 
calibration. 

2. Obtain NOAA-SR data. 

3. Apply results from event 1 to the satellite data to obtain 
daily reflected flux to space. 

4. Compute minimum albedo which is a nearly cloud-free albedo, 

5. Compute monthly average albedo. 

6. Complete development of a cloudiness parameter (effective 
could amount) which is to be derived from daily albedo and 
minimum albedo, 

7. Obtain global solar pyranometer data for selected sites in 
the continental United States, Data from 9 sites with newer 
model Eppley pyranometers were received. 

8. Compute shortwave flux absorption in the atmosphere over the 

9 pyranometer sites and normalize this absorption to a fractional 
absorption. 

9. Obtain daily RAOB data for all United States' stations. 

10. Compute daily opti-d path length due to water vapor for each 

RAOB station. 
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11. Parameterize atniospheric absorption over the nine surface 
pyranometer sites using satellite derived cloud parameter and 
optical path length results. 

12. Begin research on how the atmosphere absorption parameterization 
is to be partitioned in space. The parameterization will differ 
for each of the 9 sites. Some alternatives are to group 
parameterizations by geography* by air mass type, or combine 

all of them together into one parameterization. The last 
alternative may prove to be best with implementation of our 
new cloudiness parameter. 

13. Complete research on atmospheric absorption parameterization. 

14. Obtain a complete set of station RAOB data for the entire 
continental United States. (Same as 9 above). 

15. Compile monthly mean "normal" optical path length due to water 
vapor for each station. 

16. Obtain terrain height data for every five minutes of latitude 
and longitude. 

17. Interpolate terrain height field to same grid as NOAA-SR 
satellite data. 

18. Interpolate "normal" optical path length in three dimensions 
to NOAA-SR grid using interpolated terrain height field; and 
apply mean solar zenith factor to obtain gridded optical path 
length due to water vapor. 

19. Apply atmospheric absorption parameterization to the gridded 
optical path length field and the cloudiness parameter to obtain 
atmospheric absorption at each NOAA-SR grid point. 
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20. Combine monthly mean albedo, minimum albedo, and atmospheric 
absorption to obtain mean monthly fractional surface solar 
insolation in the NOAA-SR grid. 

21. Convert fractional insolation to solar insolation received at 

; the earth's surface. This merely involves multiplying fractional 

insolation by solar insolation received at the top of the at- 
mosphere. 

STATISTICS NETWORK 

^ Statistics besides the mean can be calculated which include consecutive 

\ days below some threshold value, the range, standard deviation, and etc. 

\ 

I These calculations cannot begin, of course, until the data processing is 

* under way as depicted in the mean network of Figure 2. 

f 

^ The flow of these events can occur simultaneously with the processing 

f • in Figure 2 by applying daily values of albedo and optical path length 

of each grid point in events 20 and 19 respectively. All beginning 

events in Figure 3 pertain to Figure 2 and are discussed in the previous 

r 

section. All other events of Figure 3 are as follows: 

21. Compute monthly mean surface insolation at each grid point. 

This event, as discussed in the previous section, may be modified 
to follow after event 24 in Figure 3. 

22. Compute daily absorbed solar energy in the atmosphere by 
applying absorption parameterizations. 

23. Compute daily fractional insolation at the earth's surface at 
each grid point. 

24. Change gridded fractional insolation to gridded insolation. 

25. Compute consecutive days of below threshold statistics. 

26. Compute variance statistics. 






56 


ALL DIGITAL VIDEO IMAGING SYSTEM FOR ATMOSPHERIC RESEARCH 

(ADVISAR) 

Appendix C 

The Advisar system is an all purpose image processing system oriented 
toward research and development in the atmospheric and earth sciences 
with particular emphasis on satellite imagery. The overi'll system is 
described in terms of its five component subsystems: 

1. Computer Hardware and Peripherals 

2. Video Display Hardware and Control 

3. Digitizing Hardware and Control 

4. System Executive Software 

5. Applications Software 

1 . Computer Subsystem 

The control system consists of a Hewlett Packard 2100-A computer 
configured with 32K memory (16 bit words), 16 I/O ports, and 2 DMA 
channels. This mini-computer is responsible for both video display 
and digitizing control along with applications software execution 
(data analysis). The peripheral hardware available to the computer 
presently consists of: 

1. One 20 megabit digital disk unit for both executive/applica- 
tion software and imagery data storage. 

2. One 45 IPS 7-track tape drive (200 DPI, 556 BPI, 800 BPI). 

3. One 45 IPS 9- track tape drive (800 BPI, 1600 BPI). 

4. One 200 cpm card reader. 

5. One 2400 baud CRT and keyboard operator terminal. 

6. One hard copy and keyboard terminal. 
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The HP 2100-A is utilized as a multi -terminal system under the 
control of an HP executive monitor. Software development and video 
display operations may occur simultaneously with very little cross 
terminal interference. All software is cataloged to the disk unit 
such that user operation consists of activating dial-up program com- 
mands. Programs are developed at either of the keyboard consoles or 
at the card reader depending on program length and complexity. Much 
effort goes into modularizing software as FORTRAN subroutines such 
that redundancy is minimized and compatibility with other systems is 
increased. The two tape drives are available for both digital imagery 
imput and output. 

2. Video Display Subsystem 

The video display subsystem presently consists of a 512 x 512 spot 
Hughes scan converter (video refresh memory), a CONRAC color monitor 
and a hand operated joystick linked to a single spot flashing cursor. 
This system was integrated and interfaced to the HP computer by Spatial 
Data Inc. The HP 2100-A is responsible for all data transfer to and 
from the scan converter; both imagery data and graphic displays can be 
output to the scan converter memory. Analog controls allow an image 
to be color enhanced during the display period; black and white enhance- 
ment is available under computer control. Another control provides a 
gradient field display at the flip of a switch. An analog planimeter is 
also available. The joystick is used to communicate coordinate requests 
as applied to the displayed image while a cursor tracks the joystick 
motion. There is an option to permanently depict a cursor location. 
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3. Digitizing Subsystem 

The digitizing capability is provided by a Spatial Data 108 
control system coupled to a Telemation (model TMC-2100) vidicon 
camera. Any portion of a photo transparency (positive or negative) 
can be digitized to a resolution of 8-bits (256 levels) under computer 
control and output to tape, disk, or scan converter. An auxiliary 
B/W monitor is utilized to portray the photographic sector under- 
going digitization. 

4. System Software 

The operating environment is under the control of the HP Real ''’ime 
Executive Monitor (RTE). This system has been minimally modified to 
provide an image processing mode convenient to R & D level users. All 
applications software is cataloged by the disk file monitor such that 

it is easily modified or edited for special applications. 

* 

5. Applications Software 

Basic applications software is available on the computer disk to 
provide; 

1. Digitizing control. 

2. Image sector tape and disk I/O. 

3. Image enhancement. 

4. Graphics display. 

5. Analog to digital conversion 

Much of the software is being redeveloped for a new, video display 
system. The future system will consist of eight 512 x 512 solid state 
video refresh memories configured with dynamic enhancement tables and 
frame sequencing capability. Individual frame resolution will be 8-bits 


or 256 levels. A computer driven cursor will replace the present 
analog cursor such that coordinate communication will contain zero 
error. Since this video refresh memory will serve as an auxiliary 
bulk Image storage device to the computer* much of the present disk 
space allocated for Image storage will be released for new, permanent 
applications software. Most of the new software will consist of 
meteorological parameter extraction techniques as applied to satellite 
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iPTICAL DATA DIGITIZER AND DISPLAY SYSTEM 

DEPARTMENT OF ATMOSPHERIC SCIENCE, COLORADO STATE UNIVERSITY 
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OD® HARDWARE COMPONENTS 

A DIGITIZING SUB- SYSTEM 

B 108 CONTROLLER AND B/W MONITOR 

C DATA COLOR EDGE ENHANCER 

D SCAN CONVERTER 

E HP 2100 -A COMPUTER 

F DIGITAL DISK 

G 9 TRACK DRIVE 

H 7 TRACK DRIVE 

I CARD READER 

J SYSTEM TELEPRINTER AND KEYBOARD 
K JOYSTICK 

L COLOR DISPLAY MONITOR 

M OPERATOR CRT AND KEYBOARD 

N SOLID STATE RAPID VIDEO REFRESH MEMORIES 
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3. Yonder Haar, T. H. and J. S. Ellis, 1975: Solar energy microclimate 
as determined from satellite observations. Proceedings of the 19th 
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63 


APPENDIX E 

Testimony before the Solar Energy Hearing 
of the 

Subcommittee on Energy Production and Supply 
and the Subcommittee on Energy Research and 
Development of the Committee on Energy and 
Natural Resources, and the Select Conmittee on 
Small Business, United States Senate 

In the chair; Senator Floyd K, Haskell 

Metals Hall, Green Center, Colorado School of Mines 
Golden, Colorado 
Wednesday, June 1, 1977 

by Professor Thomas H, Vonder Haar, Head 
Department of Atmospheric Science 
Colorado State University 
Fort Collins, CO 80523 
(303) 491-8566 
r 2944 West Dean Drive -t 
l-Ft. Collins, CO 80521 J 


Senator Haskell, ladies and gentlemen: 

Thank you very much for the invitation to present a statement 
during these hearings. My testimony is primarily that of an expert 
witness regarding the transmission of solar energy through the at- 
mosphere and the distribution and amount of solar energy reaching 
the ground. For several years my colleagues and I have carried out 
research in these areas using measurements of solar energy, clouds and 
weather data obtained from weather stations at the ground and from 
small weather satellites. Our research has been supported by state 
and federal grants, most recently from NASA and from the Colorado 
Energy Research Institute. 

Homeowners and businesses considering the use of solar energy 
need to know how much solar energy is normally available at the site 
of their home or business. Our research results demonstrate that we 
can provide this information about any location in the United States. 
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The key component In our new capability Is the use of measurements 
from the new weather satellites. 

An example of a map showing the amount of solar energy reaching the 
ground during our research test period of June, 1975, Is Included as 
an attachment to this testimony. This “big picture", a map of average 
dally solar radiation reaching the ground over the entire U.S. , shows 
that during this month some typical regional values were as follows: 

Amount Available 
Solar energy 

Re gion (calories per square centimeter per day') 


I 

t. 

i 

[ 

f 

i 

f 

i. 


NW Washington 

450 

Northern Minnesota 

450 

New Jersey 

475 

Houston-Galveston 

450 

Florida 

500 

Missouri 

550 

Northern Utah 

600 

Rocky Mountain of 
Colorado and New Mexico 

700 

Northern Central California 

750 


Some of these numbers may seem surprising; for example, that northern 
California received 50?^ more solar energy than Florida during June, 1975, 
However, this emphasizes the role of variable weather conditions on the 
solar energy supply. Using our satellite-based method (together with a 
few ground stations maintained by the National Weather Service, NOAA), we 
can determine maps of solar energy reaching the ground for enough months 
and seasons to give reliable statistics on available solar energy to all 

1 2 

1 calor1e/cm day - 0.15 BTU per square foot per hour. 
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potentfaT user's. Since the weather satellites are planned to continue 
through the 1980's and beyond* we can expect to update as needed our 
Information about solar energy reaching all areas of the U,S. 

The *^econd attachment Is a part of a weather satellite photo showing 
how very detailed views of the solar energy microclimate are possible. 

It displays the western tip of Lake Superior* Including the city of 
Duluth* Minnesota* at a time when some Ice was separated from the southwest 
shore by a narrow band of open water. Also shown Is the computer output 
map of the same small region; to Illustrate how we process the satellite 
data by computer to obtain the amount of solar energy reaching the ground. 

The very high resolution satellite views are often Important to solar 
energy users because even the average solar energy supply can change as 
much as Z0% In 20 miles near cities* mountains, coasts* etc. 

Our results thus far have demonstrated that we can use new methods 
to make available Information about solar energy supply needed by solar 
energy users. They are documented In the FINAL REPORT under NASA Grant . 
NAS5-22372 entitled: Determination of Solar Energy Microclimate of the 
United States Using Satellite Data by Thomas H. Vender Haar and James S, 
Ellis, Department of Atmospheric Science, Colorado State University, Fort 
Collins, CO 80523, April, 1977. I will provide copies of this report to 
Senator Haskell and to staff members of the various subcommittees. 

Additional copies are available on request from my office. Our report 
also Includes results for October, 1975, and explains the solar radiation 
transfer methods and computer processing techniques available at Colorado 
State University. 

What is needed to fully benefit from these preliminary research results? 
How can these results help us reach President Carter's stated goal of using 
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solar energy In 2.5 million homes by 1985? What should Congress do? In 

response to these questions posed by Senator Haskell, I conclude my 

testimony with several recommendations: 

1) We should expand the preliminary results with a renewed 

research effort to assess the detailed capabilities of using 
weather satellite data to determine the solar energy reaching 
the ground over the United States. In order to provide 
definitive information about the solar energy supply during 
all seasons; about the best processing methods; and to deter- 
tiiine the uncertainty level of the results, a good deal of 

data must be processed. I estimate the need for a two year 

effort costing $400,000 to $500,000; it could begin immediately. 
Based on our preliminary results, the probability of success 

is high and the results of the two-year study ending in late 
1979 would be of immediate use to solar energy users. A 
publication containing solar energy supply statistics and 
maps could be distributed to homeowners and businessmen through 
the government printing office, university extension offices, 
and by solar energy businessmen. 

2) Appropriate energy research groups should determine the economic 
value of updated solar energy supply data and, possibly, solar 
energy supply forecasting. 

3) Based on results from (1) and (2) above, we should consider the 
establishment of an operational solar energy supply monitoring 
program, using both satellite and ground-station measurements, 
under the auspices of the appropriate federal agency. 

End of testimony by T. H. Vonder Haar 

ATTACHMENTS (2 figures) 
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MEAN DAILY SOLAR RADIATION 
JUNE 1975 



Figure 1. Mean daily solar radiation reaching the ground over 
the Unit(5d States during June 1975 as determined from weather satellite 
observations. ^Values arc in calories per square centimeter per day. 

(1 calorie/cm day ■= 0.15 BTII per square foot per hour) 

Taken from the Final Report for MSA Grant NAS5-22372, "Deter- 
mination of Solar Energy Microclimate of the United States Using 
Satellite Data" by Thomas Vender Haar and James Ellis Tf the Department 
of Atmospheric Science, Colorado State University, Ft, Collins, CO April, 1977 
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Figure 2. Photograph of the western Lake Superior region from DMSP satellite and a digi 
tal representation of the rectangularly outlined area. 


Taken from "Solar Energy Microclimate as Determined from Satellite 
Observations" by T.H. Vonder Haar and J.S. Ellis of the Atrao'^^heric 
Science Department at Colorado State University, Fort Collins, CO; 
presented at the meeting on Optics in Solar Energy Utilization, 

August 21-22, 1975, San Diego, California. 



